Subsystem codes with spatially local generators 
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We study subsystem codes whose gauge group has local generators in the 2D geometry. It is 
shown that there exists a family of such codes defined on lattices of size L x L with the number of 
logical qubits k and the minimum distance d both proportional to L. The gauge group of these codes 
involves only two-qubit generators of type XX and ZZ coupling nearest neighbor qubits (and some 
auxiliary one-qubit generators). Our proof is not constructive as it relies on a certain version of 
the Gilbert- Varshamov bound for classical codes. Along the way we introduce and study properties 
of generalized Bacon-Shor codes which might be of independent interest. Secondly, we prove that 
any 2D subsystem [n, k, d] code with spatially local generators obeys upper bounds kd = 0(n) and 
d 2 = 0(ri). The analo gous upper bound proved recently for 2D stabilizer codes is kd 2 — 0(n). Our 
results thus demonstrate that subsystem codes can be more powerful than stabilizer codes under 
the spatial locality constraint. 



I. INTRODUCTION 

Fault-tolerant quantum information processing based 
on 2D topological quantum codes has received a con- 
siderable attention lately since it can be implemented 
on quantum machines with a geometrically local archi- 
tecture. The potential of topological codes as a viable 
alternative to concatenated quantum codes was first re- 
alized by Dennis et al [![. It was shown in [l[ that an 
active error correction in the 2D toric code permits re- 
liable storage of a logical qubit if the error rate in the 
quantum hardware is below the threshold value about 
1%. The threshold for storage of a qubit has been re- 
cently improved by Andrist et al by using topological 
color codes instead of the toric code. The success of 
topological codes was extended from a storage of a qubit 
to the universal quantum computation by making use of 
the powerful technique known as code deformations d-Q . 
These recent developments demonstrate that topological 
codes provide an attractive framework for design of new 
fault-tolerant protocols. 

In order to better understand the potential of topologi- 
cal codes for storing and manipulating quantum informa- 
tion, it is desirable to derive fundamental bounds on the 
parameters of quantum codes that stem from the spatial 
locality constraint and find families of codes that achieve 
these bounds. A progress in this direction has been re- 
cently made for 2D stabilizer codes [§-[l2j|. Such codes 
can be denned in a system of n physical qubits occupying 
sites of a regular square lattice of size yfn x y/n. Quan- 
tum information is encoded into a codespace C spanned 
by common eigenvectors of pairwise commuting n-qubit 
Pauli operators Si, ... , S m known as stabilizers, that is, 



C = {#) G (C 2 )®" : S a \r/>) = |V) for all a}. 

The locality condition is imposed by demanding that each 
stabilizer S a acts non-trivially only on a constant num- 
ber of qubits located within constant distance from each 
other. A code has k logical qubits if dim£ = 2 . It was 



shown in that any 2D stabilizer code obeys a bound 

kd 2 = 0(n), (1) 

where d is the minimum distance of a code, i.e., the min- 
imum weight of a Pauli operator commuting with all sta- 
bilizers and implementing a non-trivial transformation on 
L. The bound Eq. (p} is tight in the sense that for any 
given k and d one can construct a 2D stabilizer [n, k,d] 
code with n = 0(kd 2 ), see [ll| for details. 

The main conclusion of the present paper is that the 
bound Eq. (JTJ) can be violated in a dramatic way for 2D 
subsystem codes. Recall that a subsystem code [l3|, [l4| 
can be regarded as an ordinary stabilizer code in which 
only part of the logical qubits is used to store information. 
Accordingly, the codespace of subsystem codes can be 
decomposed as 



£ = £ 



logical 



c 



gauge 7 



dim£i ogica i — 2 , 



where £i og icai is the logical subsystem used to store quan- 
tum information, while £ ga uge represents the unused log- 
ical qubits usually called gauge qubits. The distance d 
of a subsystem code is the minimum weight of a Pauli 
operator commuting with all stabilizers and acting non- 
trivially on the logical subsystem Aogicai- The presence 
of the unused gauge qubits provides much more flexibility 
in the design of fault-tolerant gates and the error correc- 
tion for subsystem codes since one does not need to worry 
how a particular computational operation or an error af- 
fects the gauge qubits, see for instance [H| HH. By the 
same token, one should expect that spatial locality con- 
straints lead to less severe restrictions on the parameters 
of subsystem codes. 

One can characterize a subsystem code by its gauge 
group Q that includes all stabilizers and all logical opera- 
tors on the unused logical qubits, see Section [TT] for more 
details. We shall study 2D subsystem codes for which 
the gauge group has spatially local generators, that is, 
Q = (Gi, . . . ,G m ) where each generator G a acts non- 
trivially only on a constant number of qubits located 
within constant distance from each other (as for stabi- 
lizers S a , they may or may not be spatially local). For 
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such codes the tradeoff between n and d was character- 
ized in Q by showing that 



0(Vn)- 



(2) 



This bound is tight since the 2D Bacon-Shor code [l3[ 
has parameters d = yfn and k = 1. The question that 
remained open is whether 2D subsystem codes may have 
a better scaling of k compared with 2D stabilizer codes. 

In the present paper we answer this question in positive 
by proving that there exist a family of 2D subsystem 
[n, k, d) codes with both k and d proportional to y/n. 
More precisely, we prove the following. 

Theorem 1. Let a > and < f3 < 1/2 be any con- 
stants such that a+H2{(3) < 1, where i?2(/3) is the binary 
entropy. Then for all sufficiently large integers m and 
for all k < am there exists a 2D subsystem [2m 2 ,fc,<i] 
code for some d > j3m. The gauge group of this code 
has two-qubit generators of type XX and ZZ coupling 
nearest-neighbor qubits and some one-qubit generators. 

In contrast, all previously known 2D subsystem codes 
with the distance proportional to y/n, such as the 
2D Bacon-Shor code [13[ or the topological subsystem 
codes [l7( encode only k = 0(1) qubits. 

The proof of Theorem Q] relies on generalized Bacon- 
Shor codes that we introduce in Section IIII1 One can 
define a generalized Bacon-Shor code for any binary ma- 
trix A by placing physical qubits at the cells of A for 
which Ai.j = 1 and leaving the remaining cells empty. 
The gauge group Q has generators of two types: each pair 
of qubits c, d located in the same row of A contributes 
a generator X c X c i and each pair of qubits c, d located 
in the same column of A contributes a generator Z r Z c i. 
We show that the code Q has parameters [n,k,d], where 
n is the number of non-zero matrix elements in A, k is 
the binary rank of A, while d is determined by the mini- 
mum distance of the classical codes spanned by columns 
and rows of A, see Theorem [5] in Section IIIII Then we 
employ the Gilbert- Varshamov bound to prove existence 
of binary matrices with the desired properties, see Theo- 
rem [3] in Section llVI Finally, we show how to transform 
any generalized Bacon-Shor code into the spatially lo- 
cal form by introducing ancillary qubits and simulating 
each long-range generator by a chain of nearest-neighbor 
couplings, see Section IVl 

Our second result is a new upper bound on the pa- 
rameters of 2D subsystem codes whose gauge group has 
spatially local generators, namely, 



kd = 0{n). 



(3) 



It can be regarded as a generalization of Eq. (TTJ to subsys- 
tem codes. This bound is tight up to a constant factor 
since one can achieve a scaling k ~ d ~ \/n, see The- 
orem [TJ We also prove that the original bound Eq. (TTJ) 
holds for any 2D subsystem code in which both stabilizer 
group and the gauge group have spatially local genera- 
tors. The topological subsystem codes of [13, 13 provide 
an example of such codes. 



The proof of the bound Eq. ([3]) presented in Sec- 
tions IVIIVIIIVIIII requires some heavy machinery that 
builds upon techniques developed in 0, [111, [l2[ ■ Our first 
tool is the identity relating the number of logical oper- 
ators supported in two complementary regions of a lat- 
tice, see Lemma in Section fVTl It was originally proved 
by Yoshida and Chuang 12] for stabilizer codes. In the 
present paper we generalize this identity to subsystem 
codes using techniques of Our second tool is what 
we call a holographic principle for error correction, see 
Section IVIII It asserts that a non-trivial logical opera- 
tor cannot be supported in a region whose perimeter is 
smaller than the distance of the code. The analogous re- 
sult was proved in 11] for 2D stabilizer codes although 
the proof given in [Iff cannot be generalized to subsys- 
tem codes. In Section I Villi we combine these technical 
tools to prove the upper bound Eq. ([3]). 

In Section IIXI we summarize our results and discuss 
some open problems such as possible extensions of our 
constructions to 3D subsystem codes. 



II. STABILIZER AND SUBSYSTEM CODES 

The purpose of this section is to summarize the neces- 
sary facts pertaining to stabilizer and subsystem codes. 
The main idea of stabilizer codes is to encode k log- 
ical qubits into n physical qubits using a codespace 
C C (C 2 )®" spanned by states that are invariant 
under the action of a stabilizer group 5, 

£ = {|V) e (C 2 )®" : P\i)) = \if>) yPeS}. 

All stabilizers P <E S must be Pauli operators, that is, n- 
fold tensor products of the single-qubit Pauli operators 



X 



1 

1 



Y 



-i 

1 



Z 



1 

-1 



and the identity operators /. Such tensor products gen- 
erate the Pauli group on n qubits, 



V = {-yP 1 <E)P 2 



Pn}, 



where P a € {/, X, Y, Z} and 7 € {1, —1, i, —i} is an over- 
all phase factor that we shall often ignore. Non-trivial 
stabilizer codes correspond to Abelian stabilizer groups 
ScP such that -/ £ S. 

Logical operators of a stabilizer code S are Pauli op- 
erators preserving the codespace C. Equivalent ly, logical 
operators are Pauli operators commuting with every ele- 
ment of iS. Such operators generate the centralizer of S 
in the Pauli group, 

C(S) = {P eV : PQ = QP V Q e S}. 

One can always decompose the centralizer as 

C(S) — (S,Xi,Zi,..., Xk,Zk), 
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where Xi,Zi are the logical Pauli operators, while ele- 
ments of S correspond to the logical identity operator. 
Here and below we use the notation (. . .) for a subgroup 
generated by a family of operators. 

For any Pauli operator P € V let Supp(P) C A be 
the support of P, that is, the subset of qubits on which 
P acts non-trivially (by X, Y, or Z). We shall use the 
notation \P\ — |Supp(P)| for the weight of P, that is, the 
number of qubits in its support. 

The minimum distance d of a stabilizer code S is de- 
fined as the minimum weight of a non-trivial logical op- 
erator, that is, 

d = min |P|. 

pec(s)\s 

We shall use a notation [n, k, d] for a stabilizer code en- 
coding k logical qubits into n physical qubits with the 
distance d. 

To define a subsystem [n, k, d] code it is convenient to 
start from a stabilizer code S with k + g logical qubits 
for some g > 0. Let Xi, Zi, i = 1, . . . , k be the logical 
Pauli operators on the first k logical qubits that will be 
used to encode information. The remaining unused logi- 
cal operators Xi, Zi, i = k + 1 , . . . , k + g. together with 
stabilizers S generate the gauge group [14[ 

Q = (S, Xk+l, Zk+l) • • • > Xk+g, Zk+g)- 

We will assume that the n physical qubits live at vertices 
of a regular 2D lattice A of size y/n x y/n with open or 
periodic boundary conditions. Given this 2D geometry 
we demand that the gauge group Q must have spatially 
local generators, that is, Q = {G\, . . . , G m ), where the 
support of any generator G m can be bounded by a square 
box of size r x r for some constant interaction range r. 
Given the gauge group Q ', one can compute the stabilizer 
group S using the identity 

s = gnc{g), 

where C(Q) is the centralizer of Q in the Pauli group. 
Note that S may or may not have spatially local gener- 
ators. For example, any stabilizer of the 2D Bacon-Shor 
code [l3| has weight at least y/n. On the other hand, for 
the topological subsystem codes [Ij} generators of S have 
geometry of closed loops of constant size. In general, one 
can think of generators of Q as stabilizers broken into "lo- 
cal chunks" such that an eigenvalue of any stabilizer can 
be inferred by measuring eigenvalue of sufficiently many 
generators of Q . Subsystem codes with an Abelian gauge 
group are equivalent to ordinary stabilizer codes (such 
codes have no gauge qubits and thus Q = S). 

Logical operators of a subsystem codes are Pauli op- 
erators preserving the codespace C. Equivalently, logical 
operators are elements of the centralizer 

C{S) = (g,x 1 ,z 1 ,...,x k ,z k ). 

Since the encoding is defined only modulo gauge oper- 
ators, non-trivial logical operators are elements of C(S) 



that are not in Q. In the case of subsystem codes we 
shall often use the term bare logical operators which refers 
to elements of C{Q)\Q. Bare logical operators preserve 
the codespace C and act trivially on the gauge qubits. 
They should not be confused with dressed logical opera- 
tors which are elements of C(S)\Q. The identity 

c(s)=c(g)-g 

implies that any dressed logical operator can be repre- 
sented as a product of a bare logical operator and a gauge 
operator. The minimum distance of a subsystem code is 
defined as the minimum weight of a non-trivial dressed 
logical operator, 

d = min \P\ = min min IPGI. 

Pec(s)\g Pec(g)\g Gee 



III. GENERALIZED BACON-SHOR CODES 

In this section we introduce a generalization of the 2D 
Bacon-Shor code [1 3j j and describe its main properties. 
For the sake of clarity we shall ignore the issue of spatial 
locality until Section IVl 

Let A be an arbitrary matrix of size mxm with entries 
and 1. We shall label cells of A by pairs of indices 
c = Let n = \A\ be the Hamming weight of A, 

i.e., the total number of non-zero matrix elements. To 
define a subsystem code associated with A let us place 
a physical qubit at each cell c = with A^j = 1. 

Hence there are totally n physical qubits. The remaining 
cells for which Aij — are kept solely for illustrative 
purposes since there are no qubits located at these cells. 
A subsystem code associated with A has a gauge group 
Q generated according to the following rules: 

• Every pair of qubits c, c' located in the same row 
of A contributes a generator X C X C > 

• Every pair of qubits c, c' located in the same column 
of A contributes a generator Z c Z c i 

In the special case when Aij = 1 for all the code 
Q coincides with the standard 2D Bacon-Shor code [l3j . 
Consider as an example a binary matrix 



.4 : 



1 1 

1 1 

1 1 



The corresponding subsystem code has n — 6 physical 
qubits. The corresponding gauge group Q has generators 

Q = (^1,1^1,2) 2^2,3, ^3^1X3^, 
Zl^Zs,!, ^1,2^2,2) -^2, 3^3, 3)- 

Let us explain how to relate the parameters of the code Q 
to certain algebraic properties of the matrix A. It will be 
convenient to introduce a linear subspace C co \ C {0, l} m 
spanned by the columns of A and a linear subspace 
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C r ow C {0, l} m spanned by the rows of A (throughout 
this section all linear subspaces are defined over the bi- 
nary field F2). One can regard C co i and C row as classi- 
cal codes encoding k — rank (A) bits into m bits. Here 
rank (A) is the rank of A over the binary field Fa. In 
the example of the 2D Bacon-Shor code the matrix A 
has rank 1 and C co i, C row are one-dimensional subspaces 
that include only the all-zeros and the all-ones vectors. 
Let cZ co i and d row be the minimum Hamming weight of 
non-zero vectors in C co \ and C row respectively. In other 
words, d co i and e? row is the minimum distance of the clas- 
sical code C C oi and C row respectively. 

Theorem 2. Let A be an arbitrary binary matrix and 
Q be the subsystem code associated with A as described 
above. Then Q encodes k = rank (A) qubits into n = \A\ 
qubits with the minimum distance d — min (d mw , d co \). 

Before we proceed with the proof of the theorem let 
us make several remarks concerning the definition of Q . 
Firstly, the set of generators of Q introduced above might 
be overcomplete. For example, suppose c, c', c" is a triple 
of consecutive qubits that belong to the same row of A. 
Then clearly X C X C » = (X c X c r)(X c 'X c >r), that is, it suf- 
fices to retain only the generators X c X c i and X c >X c n. In 
general, it suffices to retain generators X c X c i and Z c Z c i 
for consecutive pairs of qubits c, c' that belong to the 
same row and the same column of A respectively. Sec- 
ondly, the generators of Q are not necessarily spatially 
local. For example, if some row of A contains an isolated 
l V separated by a long string of zeros on the left and on 
the right, there is no choice of spatially local generators 
for Q. We shall explain how to circumvent with diffi- 
culty in Section |V] by placing two ancillary qubits into 
each cell of A with Aij = and splitting each long-range 
generator into a chain of short-range generators. Finally, 
let us point out that a similar but technically different 
construction of subsystem codes was described by Bacon 
and Casaccino in The construction of [19| starts 

from a pair of classical linear codes C\ = [n\, k±, d\] and 
C2 = [^2,^2,^2]- A quantum subsystem code is then 
defined by placing a physical qubit at every cell of a ma- 
trix A of size n\ x n-i. The X-part of the gauge group is 
defined by replicating the parity checks of Ci in every col- 
umn of A (in the X-basis). Similarly, the Z-part of the 
gauge group is defined by replicating the parity checks 
of C2 in every row of A (in the Z-basis). The resulting 
subsystem code has parameters \n\nn, k\ki, min (di, cfc)]. 
The main difference between our construction and the 
one of is that we allow to use different classical codes 
in different rows and columns of A (although each of 
these codes is simply the repetition code on some subset 
of qubits). Also, as one can see from Theorem [21 the 
two constructions result in quantum codes with different 
parameters. 

Using Theorem [2] one can easily get upper bounds on 
the number of logical qubits k and the distance d of gen- 
eralized Bacon-Shor codes. Firstly, we claim that for any 
binary matrix A one has drow^coi < n. Indeed, since any 



column of A has weight at least d co i, the matrix A must 
contain at least d co \ non-zero rows. Each of these rows 
must have weight at least d r0 w- Hence the number of non- 
zero matrix elements in A is at least d TOV! d co \. Applying 
Theorem [2] one arrives at 

d 2 < drowrfcoi < n. (4) 

Although it is not necessary, let us point out that the 
bound Eq. ^ is not as good as it could be. In Ap- 
pendix A we shall prove a slightly stronger bound 

2d row GU(l - 2~ fc ) < n (5) 

and construct a family of binary matrices that achieves 
this bound. 

Furthermore, since A must contain at least k non-zero 
rows and at least k non-zero columns, we conclude that 
kd TOW < n and kd co \ < n. Theorem [5] then implies that 

kd < n. (6) 

In Section ITVl we shall use Gilbert- Varshamov bound to 
prove that there exists a family of codes that achieves the 
bounds Eqs. (|4I6|) up to a constant factor asymptotically 
in the limit n —> 00. The corresponding binary m x m 
matrices A are defined for all sufficiently large m and 
obey the scaling rank (A) > am, and <i co i, d row > /3m for 
some constants a, j3 > 0. It leads to the scaling k > am, 
d > j3m, and /3 2 m 2 < n < m 2 . 

Proof of Theorem [2j We shall use notations X^j and 
Zi t j for the Pauli operators acting on a qubit located at 
a cell Let us begin by describing the centralizer 

C{Q). For any row i define a row operator Ri acting by 
Z on every qubit located in the z-th row: 

Ri = \\ Z it j. 

Similarly, for any column j define a column operator Cj 
acting by X on every qubit located in the j-th column: 

Cj = Yl Xl 'i- 

i ■ Aij=l 

Proposition 1. The centralizer C{Q) is generated by the 
row and column operators, 

C{G) = (Rl, ■ ■ ■ , Rm, Ci, . . . , C m ). (7) 

Proof. Indeed, let us check that Ri £ C{Q) for any row i. 
A generator X c X r j located at the row i anti-commutes 
with Ri at both cells c and c'. Thus X C X C > commutes 
with Ri. In addition, Ri commutes with generators 
X C X C > located at rows i' ^ i since their supports do not 
overlap. It also commutes with generators Z c Z c i since 
they are both operators of Z-type. Hence Ri € C(Q). 
The same reasoning shows that Cj € C(Q). Conversely, 
let P z G C{Q) be a Pauli operator of Z-type. If P z acts 
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by Z on some qubit c, it must act by Z on every other 
qubit c' located in the same row since P z has to com- 
mute with all generators X c X c i in this row. It shows that 
P z is a product of the row operators over some subset 
of rows. Similarly, any operator P x G C(Q) of X-type is 
a product of the column operators over some subset of 
columns. It proves Eq. (J7J. □ 

Note that the supports of Ri and Cj overlap on exactly 
one qubit if Ajj = 1 and do not overlap if A^j =0. It 
follows that the matrix A controls the commutation rules 
between the row and column operators, namely, 

RiCj = (-1) A ^ CjRi (8) 

for all pairs i, j . Using Proposition[T]we can parameterize 
any JT-type operator P x £ C(Q) by a binary string x G 
{0, l} m such that 

m 

/ ,A IK; II x id . (9) 

J=l iJ:AijXj=l 

Similarly, any Z-type operator P z £ C(Q) can be param- 
eterized by a binary string z G {0, l} m such that 

m 

P Z = Y[R?= J[ Z itj . (10) 

i— 1 i,j : ZiAij —1 

The commutation rules Eq. ([5]) then imply that 
P X P Z = ^ X f Ax P z P x , z T Ax ^.1,.,, 

(11) 

Recall that the stabilizer group S of a subsystem code is 
defined as S = Qr\C{Q). From Eq. flTTJl we infer that 
P x commutes with all elements of C(G) iff x £ Ker(A). 
In this case one has P x £ C{C{Q)) = G, that is, P x £ S. 
The same argument shows that P z £ S iff z £ Ker(A T ). 
Thus stabilizers of X-type and Z-type can be identified 
with right and left zero- vectors of A respectively. Using 
the standard Gram-Schmidt orthogonalization one can 
choose k — rank (A) pairs of operators P x ,P Z £ C(G), 
a = 1, . . . , k, such that P x are linear combinations of the 
column operators, P z are linear combinations of the row 
operators, 

P X P Z = (-l) 5 ^P z P x , 

and 

C(G) = (S,P X ,...,P X ,P Z ,...,P Z ). 

It shows that P x , P z are the bare logical Pauli operators 
and the code G has k logical qubits. 

Let us now determine the distance of G- Consider some 
bare logical operator P x defined in Eq. ©. Let us an- 
alyze how one can reduce the weight of P x by multi- 
plying it with the gauge operators. Obviously, if P x 
has even weight in some row i, that is, Y]. Aj^Xj = 



(mod 2), one can completely cancel P x in this row by 
multiplying it with the generators X c X c i located in this 
row. On the other hand, if P x has odd weight in some 
row i, that is, J2j AijXj = 1 (mod 2), the best one can 

do is to reduce the weight of P x in this row down to 
1. Also it is clear that multiplying P x with gauge op- 
erators of Z-type cannot decrease its weight. It shows 
that the minimum weight of a dressed logical operator 
corresponding to P x is equal to the number of rows i for 
which ■ AijXj = 1 (mod 2). The number of such rows 
is nothing but the Hamming weight of the vector Ax. 
Note that Ax =t whenever P x is a non-trivial logical 
operator since Ax = implies that P x £ S, see above. 
We conclude that the minimum weight of X-type dressed 
operators coincides with the minimum Hamming weight 
of a vector Ax where x ^ Ker(A). Since such vectors 
span the subspace C co i, their minimum weight coincides 
with the distance d co i. Similar arguments show that the 
minimum weight of Z-type dressed operators coincides 
with d row . □ 



IV. GILBERT- VARSHAMOV BOUND FOR 
BINARY MATRICES 

Let B(m, k) be the set of all mx m binary matrices with 
the rank k over the binary field F2. Given any matrix 
A £ B(m, k) let C co \{A) C {0, l} m be the linear subspace 
spanned by the columns of A. Similarly, let C row (A) C 
{0, l} m be the linear subspace spanned by the rows of A. 
Let d co \(A) and d Ttyw (A) be the minimum distance of the 
classical code C co i(^4) and C row (^4) respectively. We shall 
use the notation H 2 (p) = — plog 2 p — (1 — p) log 2 (1 — p) 
for the binary entropy. The main result of this section is 
the following theorem. 

Theorem 3 (Gilbert- Varshamov bound). Let a > 

and < /3 < 1/2 be any constants such that a < 

1 — H2 (/?) • Then for all sufficiently large m and for all 
k < am there exists a matrix A £ B(m, k) such that 
d C o\{A) > /3m and d r0 w(A) > /3m. 

Proof. Our arguments will rely Gilbert- Varshamov 
bound for classical codes as presented in [20]. 

For any non-zero vector x £ {0, l} m let W(x) be the 
number of matrices A £ B(m,k) such that x £ C co \(A). 
We claim that W(x) does not depend on x as long as x 
is a non-zero vector. Indeed, let e\ = (10 ... 0) be the 
first basis vector of {0, l} m . Let W be the number of 
matrices A £ B(m,k) such that e\ £ C co i(^4). Any non- 
zero x £ {0, 1}™ can be written as x = Re\ for some 
R £ B{m, m). By definition of the column space we have 
x £ C co i(A) iff x = Ay for some y £ {0, 1}™. Hence 

x £ C col (A) iff ei G C col {R- x A). 

Since the map A — > RA defines a bijection from the set 
B(m, k) to itself, we conclude that W(x) = W. Since 
the rank of a matrix is invariant under transpositions, we 
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conclude that for any non-zero x £ {0, l} m the number of 
matrices A £ B(m, k) such that x £ C Iaw (A) also equals 
W. 

Define a matrix T x ,a labeled by non-zero x £ {0, 1}™ 
and A £ B(m,k) such that 

r = f 1 if x e C co i(A), 
^ | otherwise. 

For any A 6 B(m 7 k) one has X^o^- 4 = 2 fe — 1 since 
C C oi(^4) has dimension k. Also for any i / one has 
Y,AeB(m,k) ^ X ' A = W - Hence we arrive at 

(2 k -l)\B(m,k)\ = (2 m -l)W. (12) 

Let d = \(3m~\ and N be the number of matrices A £ 
B(m,k) such that C co \(A) or C IOW (A) contains a non- 
zero vector with weight smaller than d. It suffices to 
check that for the chosen d and any k < am one has 
N < \B{m, k)\. Applying the union bound one gets 

d-l , V 

N <2W^2[ m ) < 2W ■2 mH2[(d - 1)/m) < 2W2 mH2{p) . 
i=l ^ 1 ' 

Here we used the assumption f3 < 1/2. Using Eq. (fl2)) 
and a bound (2 m — 1) _1 < 2 1 ~ m we arrive at 

N <4. 2 fc - m + miT =(« • |B(m,*)|. 

We conclude that A" < |S(to, fc)| whenever 

m(a- l + flaCS)) < -2- 

Since we assumed that a < 1 — H2{(3), it holds for all 
sufficiently large m. □ 

V. BREAKING UP THE LONG-RANGE 
GENERATORS 

Let A be an arbitrary m x m binary matrix and Q be 
the subsystem [n, k,d] code associated with A. As was 
pointed out in Section IIII1 one can choose an indepen- 
dent set of generators of Q that includes operators X C X C > 
for consecutive pairs of qubits c, c' that belong to the 
same row, and operators Z c Z c i for consecutive pairs of 
qubits c, c' that belong to the same column. We shall 
consider the 2D geometry in which cells of the matrix A 
are identified with sites of a two-dimensional lattice of 
size m x m. 

Consider a pair of consecutive qubits c, c' in some row 
i and the corresponding generator X c X c i . In general the 
cells c and c' are not nearest neighbors, so there might 
be one or several empty cells in the i-th row (those with 
Ai.j — 0) between c and c'. Let these empty cells be 
ci, ... ,c p , see Fig. [TJ 

Our strategy will be to simulate the long-range gen- 
erator X c X r j by adding an ancillary qubit at every cell 



I 
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c 


c, 


c 2 




c 4 


c 5 


C 



FIG. 1. The cells c and c' represent a pair of consecutive 
qubits in some row of A that are not nearest neighbors. The 
cells a , . . . , Cp are empty. 

Ci, . . . , Cp and connecting c, c' by a chain of short-range 
generators, 

X C X C , = (X C X C1 ) ■ (X C1 X C2 ) ■ ■ ■ (X Cp X c ,). (13) 

More formally, we shall define a new subsystem code Q' 
obtained from Q by taking out the long-range generator 
X c Xd and adding all the short-range generators that ap- 
pear in the right-hand side of Eq. (p~3|) . In addition, Q' 
will contain one-qubit generators Z C1 , . . . , Z Cp on every 
added ancillary qubit. We have to verify that the codes 
Q and Q' have the same parameters k and d. It follows 
from the following lemma. 

Lemma 1. Let Q be any subsystem [n,k,d] code. Let 
q be one of the physical qubits of Q and a be an extra 
ancillary qubit. Define a new subsystem code 

Q' = {G, X q X a , Z a ) , (14) 

where all elements of Q act trivially on the ancillary 
qubit. Then Q' has parameters [n + 1, k,d]. 

Indeed, applying the lemma with q = c and a = c\ 
we obtain a new code Q' with the same parameters k 
and d, with an extra short-range generator X c X Cl , and 
an extra one-qubit generator Z Cl . Although Q' inherits 
the long-range generator X c X c i € Q, we can now replace 
this generator by X Cl X c * = (X c X Cl ) ■ (X C X C >). Note 
that X Cl X c i has a shorter length. Now we can apply 
the lemma again with q = c\ and a = C2 which breaks 
the generator X Cl X c > into a pair X Cl X C2 , X C2 X c i and an 
extra generator Z C2 . We can continue this process until 
all cells Ci , . . . , c p are occupied by ancillary qubits and 
the long-range generator X C X C : is replaced by a chain 
of short-range generators, see Eq. (|13[) . and one-qubit 
generators Z cl ,... : Z Cp . 

Similar arguments can be applied to break up all long- 
range iv-type generators. The resulting code Q" has at 
most two qubits at every cell of A. It has two-qubit 
generators X c X c r, Z c Z c i coupling only horizontal and 
vertical nearest neighbor cells respectively. In addition, 
Q" has one-qubit generators Z c and X c for some cells c. 
Lemma [T] implies that the code Q" has the same param- 
eters k and d. 

Now we can easily prove Theorem [1] Let m be any 
sufficiently large integer. Theorem [3] implies that for any 
k < am there exists a matrix A £ B(m, k) such that 
d co \{A) > [3m and d IOW (A) > [3m. By Theorem the 
corresponding subsystem code Q has parameters [n, k, d], 
where n < m 2 and d — min (ti co i, ci r0 w) — sgg The- 
orem [5] Transforming Q into a local form as explained 
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above increases the number of physical qubits due to the 
addition of ancillas. Since each cell of A contains at most 
two qubits, we can get a [2m 2 , k,d] code. (Some cells of 
A may contain only one qubit or no qubits at all. We can 
add extra unused gauge qubits to each of those cells thus 
making the total number of qubits 2m 2 .) In the rest of 
this section we prove Lemma [T] 

Proof. Consider an auxiliary subsystem code with a 
gauge group Q — (Q 7 X a ,Z a ), where all elements of Q 
act trivially on the ancillary qubit. Clearly Q is obtained 
from Q by adding one gauge qubit and thus Q has pa- 
rameters [n + 1, k, d}. Let U be the CNOT gate with the 
control qubit a and a target qubit q, see Fig. [5] One can 
easily check that U maps Q to Q' , that is, 



g' = {p' = upu\ P e §}. 



(15) 



It follows that Q' and Q have the same number of logical 



a 



U = 



q B 



FIG. 2. The conjugation by U maps generators of G to gen- 
erators of Q' since UX a lI t = X q X a and UZ a U^ = Z a . 

qubits, that is, Q' has parameters [n + l,k, d'] for some 
distance d' . Let us show that d 1 — d. Indeed, Eq. (fl5|) 
implies that P is a dressed logical operator for Q iff P' = 
UPW is a dressed logical operator for Q' . Let P be a 
dressed logical operator of Q with the minimum weight, 
that is, \P\ = d. Since X ai Z a G G, minimality of the 
weight implies that P acts trivially on the qubit a. Then 
P' = UPU' is a dressed logical operator for Q' and P' 
acts on the qubit a either by identity or by Z. More 
specifically, P' = PZ" where a — iff P acts on q by / 
or X, and a = 1 iff P acts on q by Z or Y . Since Z a E Q' 
we conclude that P is also a dressed logical operator for 
G' and thus d' < \P\ = d. 

Conversely, let P' be a dressed logical operator of G' 
with the minimum weight, that is, \P'\ = d! . Using the 
gauge operators X q X a and Z a we can cancel the action 
of P' on the qubit a without increasing its total weight, 
so we can additionally assume that P' acts trivially on 
a. Then P = U'P'U is a dressed logical operator for G 
and P = P'Z^ for some a G {0, 1}. Since Z a G G we 
conclude that P' is also a dressed logical operator for G 
and thus d < \P'\ = d! . We have shown that d' = d. □ 



VI. THE UPPER BOUND: TECHNICAL TOOLS 

Let us begin by introducing some more notations. We 
shall assume that the n physical qubits occupy sites of a 



2D lattice A of size yfn x y/n. For any subgroup G Q V 
and any subset of qubits M C A introduce a group 

G(M) = {PeG : Supp(P) C M] 

which includes all elements of G whose support is con- 
tained in M. In particular, 'P(M) is a group of all Pauli 
operators whose support is contained in M. For any sub- 
set M C A let M = A\M be the complement of M. 
Introduce also a group 

Gm = {Pe V(M) : PQeG for some Q G V(M)} 

which includes all Pauli operators P G V{M) that can 
be extended to some element of G ■ In other words Gm is 
a group obtained by restricting elements in G to M. By 
definition G(M) C Gm ^V{M)<ZV. 

Throughout this paper we shall ignore overall phase 
factors of Pauli operators. Then the Pauli group V can 
be regarded as the 2n-dimensional binary space [2l| such 
that multiplication of the Pauli operators corresponds to 
addition of the binary strings modulo two. For any sub- 
group G C V define its dimension dimC/ as the smallest 
number of Pauli operators generating G ■ In particular, 
dimP = 2n and dim V(M) = 2\M\. The standard stabi- 
lizer formalism [2l| implies that 



and 



C(C(G)) = G 



dimC(C?) + dim G = 2n 



for any subgroup G ■ We shall need the following simple 
fact. 

Fact 1. Let G Q V be any subgroup and M C A be any 
subset of qubits. Suppose 2 1 M \ < dim G- Then G contains 
at least one non-identity operator acting trivially on M . 

Proof. Indeed, the constraint that P acts trivially on M 
leads to a system 2\M\ binary linear equations which has 
a non-trivial solution whenever 2\M\ < dim (J. □ 

Consider a subsystem code with a gauge group G and 
a stabilizer group S = G CiC(Q). For any subset of qubits 
MCA let l(M) be the number of independent dressed 
logical operators supported inside M, that is, 

l(M) = dimC{S M ) n T(M) - dim G(M). 

Let Zb a re(-^) be the number of independent bare logical 
operators supported inside M, 

Zbarc(M) = dimC(G M ) n V(M) - dimS(M). 

Our first technical tool will be the following lemma which 
was originally proved for stabilizer codes by Yoshida and 
Chuang [3. 

Lemma 2. Suppose a subsystem code G has k logical 
qubits. Then 

l h& r C {M)+l(M)=2k 

for any subset of qubits MCA. 
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It can also be regarded as a generalization of the 
"Cleaning Lemma" proved for subsystem codes in 
The Cleaning Lemma of Q dealt only with a special 
case Zbare(-^0 = or l(M) — 0. Specializing Lemma [5] 
to stabilizer codes (Q = S) one gets a simpler statement 
l(M) + l(M) = 2k which coincides with the result proved 
in [12j ■ Note that Lemma [5] does not need any spatial 
locality properties. 



Lemma 3 (Restriction Lemma) . Suppose a gauge group 
Q has spatially local generators with an interaction range 
r. Choose any subset M C A and consider the subsystem 
code with the gauge group Qm ■ Then one of the following 
is true: 

(1) The code Qm has no logical qubits, 

(2) The code Qm has distance at least d— \dM\ 
where d is the distance ofQ. 



Proof. Let m = lh&re(M) and let Pi, ... , P m be m inde- 
pendent bare logical operators supported inside M, that 
is, 

P a eC{Q)\Q and Supp(P a ) C M (16) 

for all a — 1, . . . , m. Since the code has k logical qubits, 
one can choose t = 2k — m independent bare logical oper- 
ators Qi,...,Q t € C(Q)\Q commuting with Pi, . . . ,P m . 
Then we have 

C(Q M )nV(M) = (S(M),P u ...,P m ). (17) 

and 

PaQb = QbPa, Vtt,6. (18) 

Taking the centralizer of both parts of Eq. (fl7|) one ar- 
rives at 

Qm = V(M) n C(S(M)) n C(Pi) n . . . n C(P m ). (19) 

Represent Q a = Q' Q n Q° ut , where Q^ 1 and Q° ut are the re- 
strictions of Q a onto M and M respectively. Combining 
Eqs. (|16I18I19|) and taking into account that Pi, ... , P m 
have support only on M one gets 

QaeQ M for aUa = l,...,i. 

By definition of Qm, there must exist gauge operators 
G a e Q extending Q™, that is, Q' a = Q a G a € V(M) for 
all a = 1, . . . , t. It follows that Q' 1 ,...,Q' t are indepen- 
dent dressed logical operators supported on M, that is, 
l(M) > t. On the other hand, any dressed logical op- 
erator supported on M must commute with Pi, ... , P m , 
and thus l(M) < t. We have proved that l(M) =t. □ 

Our second technical tool is the "Restriction Lemma" 
proved in [H which relates the distance of a subsystem 
code Q defined on the entire lattice A to the distance of 
a code Qm obtained by restricting Q onto some subset of 
qubits M C A. Note a generators of Qm can be chosen 
as generators of Q restricted to M. Assuming that Q has 
spatially local generators with some interaction range r, 
the same is true for Qm- 

Definition 1. Given an interaction range r and a subset 
M C A let DM C M be the set of all sites in M that lie 
within distance r from M . 



VII. HOLOGRAPHIC PRINCIPLE FOR ERROR 
CORRECTION 

Let Q be the gauge group of some subsystem [n, k, d] 
code. We shall fix some choice of spatially local gen- 
erators of Q such that each generator has support in a 
box of size r x r for some interaction range r = 0(1). 
By definition of the distance d, no subset M C A of 
less than d qubits can support a dressed logical operator, 
that is, l(M) — whenever \M\ < d. Accordingly, one 
can choose R proportional to Vd such that no square box 
of size Rx R supports a dressed logical operator. In this 
section we derive a much stronger condition which can 
be regarded as an analogue of the famous holographic 
principle. 

Lemma 4. One can choose R = fl(d) such that no 
square box of size R x R supports a dressed logical op- 
erator. 

Loosely speaking, the lemma asserts that a non-trivial 
dressed logical operator cannot be supported on a region 
whose perimeter is smaller than the distance d (with some 
constant coefficient depending on r) even if the number 
of qubits in the interior of the region is much larger than 
d. A similar result was obtained in [ll[ for 2D stabilizer 
codes. We shall begin by proving an auxiliary lemma. 

Lemma 5. Let A, B be any disjoint subsets of qubits 
such that 1(A) = and \B\ + \d~A\ < d. Then l(AB) = 0. 

Proof. Let C be the complement of AB such that A = 
ABC. Applying Lemma [2] to the subset A we conclude 
that lb aYe (BC) — 2k. Hence we can choose a complete set 
of 2k bare logical operators X\,Z\,..., X^,Z^ € C(Q)\Q 
supported inside BC. Consider a subsystem code spec- 
ified by the gauge group Qbc that involves only qubits 
of BC. We claim that Qbc has k logical qubits and 
X a ,Z a are the bare logical operators of Qbc- Indeed, 
since X a , Z a have support on BC and commute with Q 
we infer that X a ,Z a £ C(Qbc) for all a. Accordingly, 
X a ,Z a ^ Qbc since X a Z a = —Z a X a for all a. It shows 
that Qbc has at least k logical qubits. Conversely, let 
P € C(Qbc)\Qbc be any bare logical operator of Qbc- 
By definition of Qbc it implies that P e C(Q) and P Q, 
that is, P must be a bare logical operator of the origi- 
nal code Q in which case it can be expressed in terms of 
X a ,Z a . We have shown that Qbc has k logical qubits 
and X a , Z a are the bare logical operators of Qbc- Apply- 
ing the Restriction Lemma to the code Q bc w e conclude 



9 



that Qbc has distance d! > d-\d(BC)\ =d-\dA\. The 
assumptions of the lemma then imply d' > \B\. Thus the 
code Qbc has no dressed logical operators supported in- 
side B, that is, l'(B) = (here and in the rest of the proof 
all quantities labeled by a prime refer to the code Qbc)- 
Applying Lemma[2]to the code Qbc an d the subset B we 
infer that ^are(C) = 2k. Hence we can choose a complete 
set of bare logical operators x' a , z' a £ C(Qbc)\Qbc sup- 
ported on C. But then X a , Z a are also bare logical oper- 
ators of the original code Q which implies ^bare(C) = 2fc 
for the code Q. Applying Lemma [2] to the code Q and 
the subset C we arrive at l(AB) = 0. □ 

Now we are ready to prove Lemma [4] 

Proof. Choose any R such that rR <C d. For any square 
box M of size Rx R consider a sequence of square boxes 
A\ C A2 C . . . C A p = M such that A\ has cardinality 
I Ai| < d and Aj+i is the smallest box that contains Aj 
and the boundary of Aj. Let B{ = Aj+i\Aj. Then |£j| + 
|<9A~| < 0(l)rR < d for all i = I,..., p. Since |A X | < d 
we have £(Ai) = 0. Applying Lemma [5] inductively with 
A = Ai and B = Bi we arrive at l(A p ) = l(M) = 0. □ 



VIII. PROOF OF THE UPPER BOUND 

Now we are ready to prove the upper bound Eq. 
By assumption, the support of any generator of the gauge 
group Q can be covered by a square block of size r x r for 
some interaction range r = 0(1). Consider a partition 
of the lattice A = AB shown on Fig. [31 The region A 
consists of square blocks Ai, . . . , A m of size R x R with 
R = Q(d) such that l(Ai) = 0, see Lemma 2] We choose 
the separation between adjacent blocks in A at least r 
such that any generator of Q overlaps with at most block 
in A. We claim that 



«bare(A) = 0. 



(20) 



Indeed, suppose P £ C(Q)\Q is a bare logical operator 
supported on A. Let Pi be the restriction of P onto a 
block Ai such that P = P1P2 • • ■ P m - Since any generator 
of Q overlaps with at most one block in A, we have Pi £ 
C(Q) for all i. However there must exist at least one block 
Ai such that Pi £ Q since otherwise P £ Q. Then for 
such a block we have Pi £ C(Q)\Q, that is, Pi is a bare 
logical operator supported inside A4. But this implies 
l{Ai) > /bare(Ai) > which is a contradiction. It proves 
Eq. (j20j) . Applying Lemma [2] we get 

1(B) = 2k. 

However, a subset B can support at most 2|£?| indepen- 
dent Pauli operators which implies 2\B\ > 1(B), that is, 
\B\ > k. Simple algebra shows that |£?| = 0(n/R) = 
0(n/d) and thus kd = 0(n). 

Let us now prove the stronger bound Eq. ([1} assum- 
ing that both Q and S have spatially local generators 



with a constant interaction range r and r s respectively. 
Consider a partition of the lattice A = ABC shown on 
Fig. 01 The regions A, B consist of blocks Ai, . . . , A m 
and B%, ... , B m respectively of size Rx R with R = Sl(d) 
such that l(Ai) — and l(Bi) = 0, see Lemma [4] The 
region C consists of disks of radius max {r, r s } so that 
adjacent blocks in A and adjacent blocks in B are sepa- 
rated from each other by distance max{r, r s }. Then we 
can choose generators in Q and S such that any generator 
overlaps with at most one block in A and with at most 
one block in B. Applying the same arguments as above 
we get Zbare(A) = and thus Lemma [5] implies 



Let us assume that 



l(BC) = 2k. 



|C| < k 



(21) 



and show that it leads to a contradiction. Indeed, 
choose any set of 2k independent dressed logical oper- 
ators Pi, ... , P-2k £ C(S)\Q supported inside BC and let 
Q = (Pi, . . . ,P2k)- Applying Fact Q] to region C and the 
group Q we conclude that there exists at least one non- 
trivial dressed logical operator P £ C(S)\Q supported 
only inside B. Let Pi be the restriction of P onto a 
block Bi such that P = P1P2 ■ • • P m - Since any genera- 
tor of S overlaps with at most one block in B we con- 
clude that Pi £ C(S). However, there must exist at least 
one block Bi such that Pi ^ Q since otherwise P £ Q. 
Then Pi is a non-trivial dressed logical operator, that is, 
l(Bi) > which is a contradiction. Hence Eq. (|2"Tj) is im- 
possible and we have |C| > k. Simple algebra shows that 
|C| = 0(n/R 2 ) = 0(n/d 2 ) which yields kd 2 = 0(n). 




FIG. 3. The partition of the lattice A = AB used to prove 
the bound kd = 0(n). The region A consists of square blocks 
Ai, . - . , A m of size Rx R with R = O(d) such that Z(A) = 0. 
Adjacent blocks are separated from each other by distance 
r = 0(1). It implies l ha , IC (A) = and thus 1(B) = 2k. This is 
possible only if \B\ > k which yields kd = 0(n). 



IX. CONCLUSIONS AND OPEN PROBLEMS 

In this paper we have studied subsystem codes for 
which the gauge group has spatially local generators in 
the 2D geometry. It was shown that the parameters 
[n, k, d] of such codes must obey an upper bound kd = 
0(n). We have also introduced a family of codes, the 
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formed into the spatially local form by introducing an- 
cillary qubits and simulating every long-range generator 
by a chain of short-range generators as described in Sec- 
tion [V] Finding the optimal scaling of d and k for such 
generalized 3D Bacon-Shor codes is an interesting open 
problem. 



FIG. 4. The partition of the lattice A = ABC used to prove 
the bound fed 2 = 0{n). The regions A, B consist of blocks 
Ai, . . . , A m and Si, ... , B m respectively of size R x R with 
R = Q(d) such that l(A t ) = and l(Bi) = 0. The region C 
consists of disks of radius max{r,r s } so that adjacent blocks 
in A and adjacent blocks in B are separated from each other 
by distance max {r,r s }. It implies Iharc(A) = and thus 
l(BC) — 2k. If fe violates the upper bound one would have 
1(B) > l(BC)-2\C\ > 0. Assuming that 5 has spatially local 
generators this is possible only if l(Bi) > for some block Bi 
which is a contradiction. 
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APPENDIX A 



generalized Bacon-Shor codes, that achieves this bound 
with both k and d proportional to y/n. The gauge group 
of the generalized Bacon-Shor codes involves only two- 
qubit generators of type XX and ZZ coupling nearest- 
neighbor qubits (and some one-qubit generators). It fol- 
lows that the syndrome measurement for these codes re- 
quires only eigenvalue measurements for operators XX 
and ZZ on nearest-neighbor qubits. Our proof of ex- 
istence presented in Sections IIIIIIVI is not constructive 
since it requires binary matrices achieving the Gilbert- 
Varshamov bound stated in Theorem |3l On the other 
hand, one can easily show that a random m x m bi- 
nary matrix A with a fixed rank k achieves the Gilbcrt- 
Varshamov bound of Theorem [3] with probability ap- 
proaching one in the limit m —¥ oo. Therefore for finite 
sufficiently small lattice sizes one can simply choose the 
desired binary matrix A randomly (with a fixed rank), 
compute the minimum distances d co i, d row and check 
whether the Gilbert- Varshamov bound is satisfied. 

A serious drawback of the standard 2D Bacon-Shor 
code [HI that precludes it from being used in the topo- 
logical quantum computation schemes is the lack of a 
constant error threshold in the limit of large lattice 
size [H, We expect that the same drawback is 

shared by the generalized Bacon-Shor codes introduced 
in the paper, ft is therefore an interesting open problem 
whether it is possible to construct 2D subsystem codes 
with both k and d proportional to y/n which would have 
a good behavior under random uncorrelated errors. 

Finally, let us point out that our construction of the 
generalized Bacon-Shor codes naturally extends to 3D 
Bacon-Shor codes [1 31 ) . In the 3D case the binary matrix 
A should be replaced by a three-dimensional binary array 
with qubits occupying cells with = !• The corre- 

sponding gauge group Q is generated by operators XX, 
YY , and ZZ coupling pairs of qubits that differ only in 
x, y, and z-coordinate respectively. For any choice of 
the array A the resulting subsystem code can be trans- 



Thc purpose of this section is to prove the upper bound 
Eq. (0. We shall also construct a family of codes that 
achieves this bound. 

We will show that a tuple [n,k,d mw ,d co i] can be re- 
alized by some binary matrix A only if the following 
quadratic optimization problem has feasible solutions: 



r x >0 
c x > 

^ ^ f*x — ^row 
X : xy — 1 

^ ^ C% ^ <^col 
x : x-y—1 

E " 



' x^y 



\tx e s fe (22) 
Vx G S fe (23) 
Vy G E fc \0, (24) 

Vy e E fc \0, (25) 

(26) 



x,y : x-y—1 



Here r x and c x are integer- valued variables labeled by 
binary strings x € S fc = {0, l} fc . We used the notation 
x • y = y^., xtyi (mod 2) for the binary inner product. 
Hence we can get a lower bound on n by minimizing the 
quadratic function of r x , c y defined in Eq. (|26l) subject to 
constraints Eqs. (|22H25I) . 

Let us begin by deriving analogous optimization prob- 
lem corresponding to ordinary classical codes. Let G be 
the generating matrix of some classical [n, k, d] code, such 
that G has size k x n and the rows of G form the basis of 
the codespace. For any binary string x — [x±, . . . , Xk] G 
E fe let n x be the number of columns [x%, . . . , Xk) T in the 
matrix G. Any codeword can be represented as yG for 
some binary string y G S fe . One can easily check that 
the Hamming weight of yG can be expressed as 



\VG\ 



E 

X : X -y—l 



where the summation is over binary strings x G 
Hence a tuple [n, k, d] can be realized by some code iff the 
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following optimization problem has feasible solutions: 

n x > Vx e S fe (27) 
5^ n, > d Vy e S fe \0, (28) 

x : x-y—1 

n x = n. (29) 

Here we treat n x as integer-valued variables. Indeed, any 
solution {n x } can be transformed into a generating ma- 
trix G (defined uniquely up to permutation of columns) 
of size k x n. Then Eq. (|2"9"|) implies that G represents a 
[n, k, d] classical code. 

Now consider an arbitrary binary matrix A of rank 
k. Without loss of generality the first k rows and the 
first k columns of A are linearly independent (other- 
wise permute rows or columns). Let G row be the gen- 
erating matrix of a classical code spanned by the first 
k rows of A. Similarly, let G co i be the generating ma- 
trix of a classical code spanned by the first k rows of 
A T . (Note that G row and G co i may have different length 
if A is not a square matrix.) Let c x be the number of 
columns x — [x%, . . . ,Xk] T in G co i. Let r x be the num- 
ber of columns x = [xi, . . . , Xk] T in G row . The variables 
c x ,r x must obey inequalities analogous to Eq. (|27I28[) 
since we assumed that G co i and G row have distance d co \ 
and d row respectively. It yields Eqs. (|22j|25[) . It remains 
to derive Eq. (f2"rJ|). Consider any row of A that starts 
with x = [x\, . . . , It can be represented as zG row 
for some z — z(x) € S fc since by assumption any row 
of A is a linear combination of the first k rows. More- 
over, the function z(x) must be linear and invertible, that 
is, z(x) = xM for some k x k invertible matrix M. As 
before, the Hamming weight of zG row can be expressed 



as 

|^G row | ^ Tyj 

v ■ y-z=l 

where the summation is over binary strings i/g E'. Since 
the number of rows in A that start from x = [x\, . . . , Xk] 
is equal to c x we arrive at 

n =\A\= 2J Cx \z(x)G row \ = ^ c x r v 

iES 1 x,y : xM-y—1 

Here xM ■ y is the inner product between binary strings 
xM and y. Since Eqs. (|22II25I) are invariant under a 
change of variables c x — > c x m for any invertible matrix 
M, we get Eq. ([26]) . 

Now we can easily prove Eq. Let r — ^2 x ^ Q r x . 
Adding up Eq. (|24|) for all y ^ we count each r x exactly 
2 k ~ 1 times, that is, we get 

r > d ro w(2 fc - l)2 1 - fc = d row (2 - 2 1 - k ). 

Then combining Eqs. (|25I26[) we get n > d co ir which is 
equivalent to Eq. ((5|). 

Let us show that Eq. is tight. Choose any integer 
k > 1 and define a matrix A of size (2 fc — 1) X (2 fc — 1) as 
a binary version of the Hadamard matrix, 

A x ^y — x • y 

with x,y € S fc \0. Obviously, A is a symmetric ma- 
trix. Its row-space and its column-space coincide with 
the Hamming code [2 k — 1, fe, 2 ]. In particular, any 
row and any column of A have weight 2 k ~ 1 . Thus we get 
d mw — d co \ = 2 k ~ 1 and n — (2 k — l)2 fe_1 . It achieves the 
bound Eq. 
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